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Abstract — This paper studies the interference channel with 
two transmitters and two receivers in the presence of a MIMO 
relay in the low transmit power regime. A communication 
scheme combining block Markov encoding, beamforming, and 
Willems' backward decoding is used. With this scheme, we get an 
interference channel with channel gains dependent on the signal 
power. A power allocation for this scheme is proposed, and the 
achievable rate region with this power allocation is given. We 
show that, at low transmit powers, with equal power constraints 
at the relay and the transmitters, the interference channel with a 
MIMO relay achieves a sum rate that is linear in the power. This 
sum rate is determined by the channel setup. We also show that in 
the presence of abundant power at the relay, the transmit strategy 
is significantly simplified, and the MAC from the transmitters to 
the relay forms the bottle neck of the system from the sum rate 
point of view. 

I. Introduction 

The capacity of the interference channel (IC) is a thirty 
years old problem in network information theory, that is of 
practical importance as well. When more than one transmitter 
and receiver want to communicate simultaneously, interference 
Umits their communication. The rate region for the simplest 
case of two transmitters and two receivers has been thoroughly 
studied, but the problem remains open for the general case. 

Recently, some good achievements have been made in 
characterizing the degrees of freedom and achievable rate 
regions of interference networks. It was shown in J^, that 
by using a simple Han-Kobayashi scheme [2], the capacity of 
a two user interference channel can be achieved to within one 
bit. For the general case of a K-usst interference network, it 
was shown in [3] that the degrees of freedom is given by K/2, 
i.e. the capacity can be well characterized by 

|-log(l + SNR) + o(SNR), 

where the second term decreases for increasing SNR. From 
a practical point of view, it is always interesting to analyze 
the performance of suboptimal schemes. For instance, in 
the rate region of a K-nsei interference channel is analyzed 
for the case in which the interference is treated as noise. 
The optimality of treating interference as noise for the two- 
user interference channel has been analyzed in [5] , [6] , [7] , [8] . 
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Power allocation strategies for the same system have been 
analyzed in [9]. Game-theoretic aspects have been considered 
in [10]. 

Another direction in the study of the IC is the interference 
relay channel (IRC), where a relay is used to support the 
communication between transmitters and receivers. This has 
gained research interest since [11 . Recently, a communication 
scheme that achieves full degrees of freedom at high SNR was 
proposed in T2] for the interference channel with a MIMO 
relay (IMRC). In this scheme, the transmitters communicate 
with the relay in a MAC phase, then the relay broadcasts the 
received data to the receivers. This is of practical interest, since 
in practice, the relay does not have knowledge of the transmit 
signals. 

In this paper, we consider the IMRC with the communica- 
tion scheme proposed in [12J. Namely, this scheme uses super- 
position block Markov encoding, beamforming, and Willems' 
backward decoding. In spite of its complexity, this scheme 
transforms the IMRC to an IC, with channel gains dependent 
on the signal power, which simplifies the study of the IMRC. 
In 12J, some power allocation strategies are considered, but 
these power allocations are not optimal; they are of interest 
for high transmit power P, where they were used to state the 
degrees of freedom of the system. We extend the study to the 
low P case, where we study the performance of this scheme, 
and propose an (approximately) optimal power allocation. 

We give the model of the IMRC in section |ll] and describe 
the communication scheme in section |III1 Then we study its 
performance at low P in section |lVl A numerical example is 
included in section [V] Finally, we conclude with section |VI] 

II. System Model 

Figure ([T]) shows a model of the IMRC. Each transmitter 
needs to communicate with its respective receiver, and the 
relay tries to support this communication. We assume that the 
transmitters and receivers are equipped with one antenna each, 
and the relay is equipped with 2 antennas. 

We denote by xi, X2, and x/j the transmitted signals of 
transmitter 1, 2 and the relay respectively, and by yi and y2 the 
received signals at receivers 1 and 2, respectively. We consider 
zero mean, unit variance, additive white Gaussian noises at the 
receivers and the relay denoted as zi, Z2, and zr. So we can 




Fig. 1. A model for the interference relay channel 



write the input-output relations as: 

yi = hiixi + h2iX2 + h.iiiXR + zi, 

V2 = hi2Xi+h22X2+'^R2^R + Z2, 
YB. = SlRXl + g2B.X2 + Zfl, 

where for i,j G {1,2}, i ^ j, ha denotes the direct channel 
gain from transmitter i to receiver i, hij the cross channel gain 
from transmitter i to receiver j, — [gn gaY' the channel 
gain from transmitter i to the relay, and h-m = [hjn^i /ifli,2]^ 
the channel gain from the relay to receiver i. The transmitters 
have a power constraint P, and the relay has a power constraint 
Pr. We assume that the relay operates in full duplex mode, 
and has global channel knowledge. 

III. Coding Scheme 

The coding strategy considered is the one proposed in |12| . 
and we will briefly explain it in this section. We consider 
transmission over a period of B blocks, where the sources 
and the relay send sequences of S — 1 messages. If a rate pair 
(i?i,i?2) is achievable in a block, then this scheme achieves 
a rate pair -^^j i?2-^g^), that approaches (i?i,i?2) as 
_B — > cx). This coding strategy at the transmitters and the relay 
is sketched in Table Ufor the general case, and in the following, 
we explain it in more details. 

A. Encoding at the Sources 

We use super-position block Markov encoding at the sources 
[H], i.e. 



xi{b) 

X2{b) 



ui{b)+u[{b), 
U2{b) +u'2{b), 



where for user i, i e {1,2}, Ui{b) is the codeword of the 

message of block b, with power pi, and = \J^Ui{b—l), 
is the codeword of the message of the previous block b — 
1, with power p[, such that pi e [0, P[, and Pi + Pi = P- 
The transmitters use predefined messages (pi and 02 as the 
messages of block 0, i.e. mi(0) and U2(0). 



B. Decoding and Re-encoding at the Relay 

The relay uses the SDMA scheme described in 
[T3l Section 10.1]. Assuming that the decoding of messages 
ui{b — 1) and 1*2(6 — 1) was successful, the relay can subtract 
them from the received signal, and then decode the messages 
uiib) and 7/2(6) using successive interference cancellation, 
achieving rate constraints given by 
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where G = [gi^^ g2R], Kp = diag{pi,p2), and I2 is the 2 x 2 
identity matrix. 

After decoding, the relay uses multimode beamforming to 
transmit to the receivers, i.e. the relay constructs the signal 

Xi?,(^) = u'jii{b)ti + U'fl2(^)t2, 

where ti and t2 are unitary 2x1 beamforming vectors. In 
our approach, ti and t2 are chosen such that they reduce 
interference at the receivers. Let pi,p2 £ [0, 1] be the power 
trade-off coefficients at the relay, i.e. the relay splits its power 
to PiPr and P2PR for u^i(6) and u'j^2{b) respectively, such 
that pi+ p2 = 1- So 



= J^"K&),for*e{l,2}. 
y Pi 

C. Decoding at the destinations 
The received signal at receiver i for block b can be written 

as 



y,(6) = huu,{b) + {hu + J ^hl^uXib) 

y Pi 

+ hj,u,{b) + (h,, + /^h^,tj)u -(6) + z„ 

V 

with i ^ j, i,j S {1, 2}. In order to reduce interference, the 
relay chooses the beamforming vectors ti and t2 such that 



h2i 



IP2Pr,T . _ n. 

— 7— hfl,it2 - 0, 
P2 



(4) 



h,2 + J^hl,t,^0. 

Pi 



Let us denote by tio and t2o the vectors ^/ and 

nd t; 

PiPr 



^y-S-t2 respectively. Since ti and ±2 are unitary, it follows 



tiol 



P[ ' 
P2Pr 



(5) 



P2 

With Q and (|5]), we get a system of two equations with two 
unknowns for each of the beamforming vectors. Notice that 
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B-1 
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Xl 


(9^1, "i(l)) 


K(l),«i(2)) 


(ni(2),«i(3)) 




(ui{B-2),ui{B^l)) 


ui{B-l) 




(9i2,W2(l)) 


(«2(1),«2(2)) 


(m2(2),«2(3)) 




{U2(B -2),U2(B -1)) 


U2{B-1) 






(«l(l),«2(l)) 


(ni(2),«2(2)) 




{ui{B-2),U2{B-2)) 


(«i(S-l),«2(B-l)) 



TABLE I 

Sketch of the superposition block Markov coding scheme, here, (pi and <f>2 are arbitrary initialization messages known by the 

TRANSMITTERS AND THE RELAY, AND {x, y) MEANS A SUPERPOSITION OF X AND y. 



these equations do not have a unique solution. Equation (j4| 
tells us that the components of t^o are linear with respect to 
each Other, while (jS]) tells us that the beamforming vector lies 
on a circle, leading to two solutions. Solving for tio and t2o, 
we get for i ^ j, i,j e {1, 2}, 



where 



Tio = nihBj.2]J ^hfj + llhfljll^ p'^^ - hijhRj,i, (7) 

with Hi e { — 1, 1}. This gives unitary ti and t2, and satisfies 
Q. This choice of t^o reduces interference seen by the 
receivers, so then we can express yi{b) as 

2/i(6) = hiiUi{b) + {ha + h]Jjtio)w-(^) + hjiUj{b) + Zi. 

Now, the receivers can use Willems' backward decoding [14] 
to decode their signals. Starting from block B, receivers 1 and 
2 have interference free signals and can decode ui{B— 1) and 
U2{B — 1) respectively. Then, in each block b, the receivers 
subtract the akeady known signals ui{b) and U2{b) from their 
received signals before attempting to decode ui{b — 1) and 
7/2(6 — 1). Now we can express yi{b) as 



Vzib) = {ha + h%^tio)ui{b) + hjiUj{b) + z. 



(8) 



As a result, the interference relay channel transforms into an 
IC. To simplify the notation, we will use /n, /12, /21, and /22 
to denote the new channel coefficients: 



fa — ha + h^jtio, /i 



II] — hij. 



(9) 



Now we can write the obtained IC input-output equations (|8j) 
as 

yi{b) = fiiu'.ib) + fjiUj{b) + Zi, 

where fa and fji depend on the channel coefficients, pi, P, 
Pr and Pi. 

IV. Performance at low transmit power P 

We aim in this section to analyze the performance of the 
given scheme at low transmit power P. Denote the optimal 
power allocation at the transmitters for a fixed power allocation 
Pi as pi and p2, and denote the rate region achieved by 
this power allocation as TZp. Then we have the following 
proposition. 



Proposition 1: The rate region TZ of the IMRC with the 
considered scheme, at low P is given by 

n = ch\ [J Tip 

vpe[o,i] 



(6) where ch{S) denotes the convex hull of S. 



A. Treating interference as noise 

Let us assume for the moment being, that we fix a choice 
of tio and t2o, and we consider a fixed power allocation at the 
relay, i.e. fixed rii and pi. Since pi < P, we can approximate 
/ii and /22 as linear functions of pi and p2 respectively as 
follows (see details in appendix lAt 



fi-i 



El 



(10) 



where we drop the arguments of fl^\ and f^^'^ for readability. 
This approximation is needed for solving our optimization 
problem, due to the fact that the argument of the square root 
in ([7]) is not concave in p,, and hence can not be optimized 
using standard convex optimization tools (e.g. I15j ). 

The receivers in the obtained IC treat interference as noise, 
resulting in rates bounded by 



Ri < log (^1 

i?2 <1ok( 1 



\\f22\\\P-p2) 

1 + 11/12 II V 



R. 



IC 



(11) 
(12) 



B. Power allocation at low P for sum rate maximization 

Up to this point, the expressions are not low-P-specific. 
From this point on, we restrict ourself to low P. We still 
consider fixed rii and pi. Let us write the rate region for this 
scenario as 

i?i < min(i?f^^,i?(^), 
i?2 < min(i?f^^,i?f ), 

Dip ^ pMAC 

It is required to find powers pi that maximize this region. In 
the following proposition, we will specify this rate region at 
low P for fixed arbitrary pi and rii, the proof is shown in 
Appendix |B] 

Proposition 2: The rate region of the IMRC, with the 
coding scheme described in section |III] with fixed Ui and pi 



can be approximated at low P as 



Ri 
R2 



< 



< 



\Sir\\ Pi 
ln2 '■ 

ln2 ' 



(13) 



where 



Pi 



P2 = 



A2 + VAi+8HM22|PSR(M22^2*2) p 
45R(/i22;^2*2) 



Ai =25R(Aiiii^ri)-|lMiill'-||gifl.ll'.and X2 ^ 2^{p22i^l2)- 

ll/^22f -||g2K||^ 

Notice that the rate bounds in proposition |2] are linear in pi 
and p2, which are functions of ni and 712, so we have the 
following corollary. 

Corollary 1: The rate region in proposition |2] is maximized 
for a fixed arbitrary pi by choosing powers 



Pi = 

P2 = 



max pi , 

niG{-l,l} 

max p2 ■ 

n2G{-l,l} 



Plugging these powers in (fT3|) . we get the region TZp. 

C. Special Case: Pji 3> P 

In this subsection, we introduce a special case, which has 
the advantage of significantly simplifying the transmit strategy. 
Namely, we consider the case of abundant power at the relay, 
i.e. Pr ^ P. In this case, we can approximate t^o as 



l|hfl,l| 



PiPr 
P-Pi 

PiPr 
P-Pi 



where rii E { — 1,1}. It follows that the coefficients of the IC 
become 




Substituting in (pT]) and (|12l) . we get the following for R^'^ 



and Ri,^: 



IC 



Ri 



R'^ 



log I 1 
log ( 1 



dct"(H)piFfl 

|hH2|P(l + ||/2l|Pp2) 

det'(H)p2Pfl 

ihm|P(l+||/l2|Ppi) 



Rt 



AP 



If Pr is high enough, then the rates with abundant relay power 
Rf^ and i?^^ are greater than the rates at the MAC side of 



If'^^ and 

Consequently, the sum rate is determined by the MAC side 



the IMRC i?f and R2 respectively for all pi and p2. 



of the IMRC, i.e. by i?f and R^^^^, and the optimal 
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Fig. 2. Components of the beamforming vector tio plotted as a function of 
PI for P = 0.1. Tlie plot also shows our approximation for the components 
of tio. 
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Fig. 3. Plots of the functions rf'^^' and r^'-^ , showing their intersection at 
the power that delivers optimal Ri . It also shows plots of the exact expressions 
f>MAC j^W (-f-Qj arbitrary small p2 < P). 



power allocation for maximizing the sum rate in this case is 

Pi = P2 = P- 

Remark 1: The expressions in section |lll] are defined for 
Pi e [0, P[, however, they can be easily modified to include 

P^^P■ 

As a result, at high Pr, the transmitters do not need to 
use super-position block Markov encoding. Each transmitter 
sends Ui{b) in block b, the relay decodes Ui{b), and then 
sends them delayed at the next block 6+1 while still using 
multimodal beamforming. In this case, we achieve Rf^^'-^ = 
log(l + ||g,«|pP). 

V. Numerical Example 
Consider the channel with parameters 

hii = /122 = 1-2 , hi2 = /121 = 0.5, 
gifl = [0.6 1.2]^ , g2fl. = [1 0.5]^, 
hfli = [0.5 1]^ , hR2 = [1 2f, 



and assume Pr = P = 0.1. For equal power split at the 
relay, i.e. p = 0.5, the components of beamforming vector 
ti, and its approximation are shown in figure (|2|. Figure ^ 
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Fig. 4. Exhaustive searcli and approximate solution of tiie power allocation 
problem for pi (normalized to P). 
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Fig. 5. Normalized sum rates for different power allocation strategies. 



P, then we obtain a sum rate that is Hnear in P at low P. 
It follows that the normalized sum rate is a constant at low 
P, given by the channel parameters and the power split at the 
relay. 

Using super-position block Markov encoding at the sources, 
beamforming at the relay, and Willems' backward decoding at 
the receivers, the IMRC transforms into an IC. We have given 
the channel gains of this IC, as functions of the parameters of 
the system, including the powers. 

Of practical interest is the case where the relay power is 
much greater than the transmit power In this case, we have 
shown that the encoding at the transmitters becomes simpler, 
since there is no need to perform super-position block Markov 
encoding. Furthermore, the MAC from the transmitters to the 
relay forms the bottle neck for the system from the sum rate 
point of view in this case. 

Given the obtained IC, the question of the optimality of 
treating interference as noise at the receivers arises. It would 
be interesting to find conditions on this channel that allow us 
to optimally treat interference as noise. This work can also be 
extended to the high power regime, where an optimal power 
allocation that maximizes the sum rate at high transmit power 
needs to be found. 



Appendix A 
Approximations for pi < P 



shows a plot of ^ and r^. Notice that in this example, if 
we choose pi ~ 0.0583, then we maximize min(r*^"^'-^, rl*-^), 
allowing us to achieve maximum 

Figure Q shows the exhaustive search (numerical) solution 
of the power allocation problem (pi) and the approximate 
solution, both normalized to P. Finally, in figure (|5|), we 
show the sum rates (normalized to log(l + for four 

different power allocations: 

• Optimal power allocation for maximum sum rate (exhaus- 
tive search), 

• Approximate power allocation as in corollary [1] 

• Equal power allocation with pi = P2 = ^'/2, and 

• Equal power allocation with pi = P2 = ^J~P for P > 1. 
Notice that at low P, our approximation (dashed line) is close 
to the maximal sum rate, and that it is constant in that region. 
Notice also that the maximum sum rate approaches one for 
large P. The power allocation p\= P2 — VP and p\= P2~ 
P/2 give a normalized sum rate approaching zero and one 
respectively at high P which confirms results in [T2] . 

VI. Conclusion 

As a result of this work, we have obtained an approximation 
for the optimal power allocation, that maximizes the sum rate 
for the given scheme. If we consider the special case of Pr = 



Since pi < P, we can approximate in (|6]) as 

-p (l + ^) using Taylor series to the first order Moreover, 
using Taylor series, the square root term in Q can be also 
approximated as 




Remark 2: Note that this approximation is precise only 
when Pi <C P, in our case, we only know that pi < P, so this 
is a rough approximation. 

After substituting in (|6| and ([9]), we get the following expres- 
sions for /ii and f22 

hi ~ lJ'iiini,pi) + iyii{ni,pi)^, 

P2 

/22 ~ M22("2,P2) + J^22(»^2,P2)-5-, 



'Hi 



where 

^J'llinl,Pl) = /ill 

i^ii{ni,pi) = ni- 

V22{n2,P2) = -n2- 

with H = [h_Ri h_R,2] 
^je {1,2}. 



h 



fll,2'll2 



As a result of ([T4| and (pT|) we can write the rate region at 
low P as 



hR2,2 

det(H) 



+^2 



pi_PRdet(H) 
2PS'i 

hR2,2h21 
hRl,2 

det(H) 



Si - ni 



huh 



R2,l 



T'R2,2 



/92Pi?det(H) 



n2- 



ft.2l/l 



fil,l 



hRl,2 



2PS2 



S^ 



P^PR I 
P I 



Appendix B 

Low P APPROXIMATIONS 



In the following, we state the proof of Proposition |2] We 
consider low P, i.e. P ^ 0, and since pi < P, it follows that 
Pi ^ 0, i e {1, 2}. Equations ([T]) and (j2]) can be respectively 
approximated at low P as 



MAC 



llglflll Pi 

ln(2) 

l|g2flH^P2 

ln(2) 

Equation ((3]) can be re-written as 



MAC 



(14) 



MAC 



P 



MAC 



log{apiP2 + Ppi + 7P2 + 1), 



(15) 



where 



l.9ll|l 11.9221 



1.92111 113121 



/3-||gii||" + ||5i2| 

2 



7 



1521 1 



522IP = ||g2fl|P, 



and this can be approximated at low P as 



R 



MAC 



llglflllPl 



llg2fl|l P2 



ln(2) 



(16) 



Notice that the bound R 



MAC 



is redundant and needs not to 



be considered for low P. Now, equations (fTTT) and (fT2|) can 
be approximated as 



K2 



(||Aiii||2p+(25R(Ain^^i*i)-||/xiif)pi 



ln(2) 

-25R(A.iii.ri)p?/P) - 

-(llA*22||'P+(23fi(Ai22i 



M22|P)P2 



i?l 
R2 



< 
< 



MAC JC 



min(r2^^-^^,r2^^) 



In order to maximize this rate region, we would like to 
choose a power allocation that maximizes min(rf^"^'^, rC^) 
and niin(r2^"^'^, r-j*^) over pi and p2 respectively. Since 



rf"^ —^'('^ is a quadratic function of pi, and rf^'^'^—r^*^ < 



for pi = 0, r{' 



MAC 



^IC 



> for pi — P, then rf 



MAC_^IC 
1 

MAC _ ,JC 
1 



MAC 



r. 



IC 







admits a solution pi e [0, P]. Similarly, 
admits a solution p2 E [0, P]. After solving the resulting 
quadratic equations, we get 



Pi = 



P2 



4jR(Mii^ri) 
A2 + v/Ai+8||rfjR(M22^^2* ' 

45R(/i22;^2*2) 



p, 



-p 



-5l252l5ll522-5ll5225l2521: 



with Ai = 23fi(/iiiiy*i) - ||^ii||2 - ||gii^||2 and A2 = 
25R(/i22i^22) ^ 1 1/^22 IP — \\S2R.\\'^- Substituting these powers in 
(fT4l) gives us the rate region achievable by this scheme at low 
P. 
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